Anomalous diffusion, nonlinear fractional Fokker-Planck equation and solutions 
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We obtain new exact classes of solutions for the nonlinear fractional Fokker-Planck-like equation 
dtp = d x {D(x)d%~ 1 p" — F(x)p^ by considering a diffusion coefficient D = X>|x|~ e [9 G 1Z and 
V > 0) and a drift force F — —k\x + fc 7 s|a^j 7_1 (fei,fe 7 ,7 £ TV). Connection with nonextensive 
statistical mechanics based on Tsallis entropy is also discussed. 

PACS numbers: PACS numbers: 82.20.Db,66.10.Cb,05.60.+w, 05.40.+j 
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I. INTRODUCTION 



The nonlinear or fractional Fokker-Planck-like equa- 
tions [0, ||, ||, ||, ||] have been used to analyze several 
physical situations that present anomalous diffusion. A 
typical nonlinear diffusion equation is dtp — T>d 2 x p u , 
usually denominated as porous medium equation. It has 
been employed in the analysis of percolation of gases 
through porous media [y > 2) JfJ, thin saturated re- 
gions in porous media (y = 2) M, a standard solid-on- 
solid model for surface growth (y = 3), thin liquid films 
spreading under gravity (v = 4) among others ||. 
In the presence of an external drift force, a drift term is 
incorporated to the porous medium equation. For exam- 
ple, such situation has been considered in the study of 
escape time, or mean first passage time, by considering 
a nonlinear Fokker-Planck-like equation, leading eventu- 
ally to a generalization of the Arrhcnius law On 
the other hand, the fractional diffusion equations are 
frequently considered as alternative approaches to con- 
tinuous time random walk models, generalized Langevin 
equations, or generalized master equations. For instance, 
we can mention that a generalized master equation, re- 
lated to a fractional Fokker-Planck equation, has been 
used in the modelling of nonmarkovian dynamical pro- 
cesses in protein folding jll| , describe relaxation to equi- 
librium in system (such as polymers chains and mem- 
branes) with long temporal memory p2| , and in anoma- 
lous transport in disordered systems p3[. Another class 
of anomalous diffusion can be obtained using a spatial 
dependent diffusion coefficient D = D{x). For instance, 
in |fI3 D oc \x\~ 6 has been employed to analyse the dif- 
fusion in a fractal medium. 

Typical physical situations where these kinds of equa- 
tions are applied concern anomalous diffusion of the cor- 
related type (both sub- and superdiffusion; see jL5j and 
references therein) and of the Levy type (superdiffusion; 
sec jf(| and references therein). The anomalous cor- 
related diffusion, when it has a finite second moment 
(x 2 ) cx t a , corresponds to superdiffusion, normal diffu- 
sion and subdiffusion for a > 1, a = 1 and < a < 1, 
respectively; a = basically indicates a localization. For 
Levy diffusion ( Levy distributions), there is no finite 



second moment, i.e., (x 2 ) diverges. 

We dedicate the present work to investigate the fea- 
tures that emerge from considering a mixing between 
fractional derivative, nonlinear aspects and spatial de- 
pendent diffusion coefficient. More precisely, we focus 
the nonlinear fractional Fokker-Planck-like equation: 



d . . d 
M pix > t] '= d~x 



D(x)^-^[p(x,tW - F(x)p( Xl t) 



(1) 



where e K, D(x) oc \x\~® is a (dimensionless) diffu- 
sion coefficient (9 £ TV), and F(x) = —dV(x)/dx is a (di- 
mensionless) external force {drift) associated with the po- 
tential V(x). We use the Riemann-Liouville operator [|l8| 
for the fractional derivative, and we also work with the 
positive x axis and, later on, we employ symmetry to ex- 
tend the results to the entire real axis (in other words, 
we are working with In addition, the ini- 

tial condition p(x, 0) = S(x) and the boundary condition 
p(±oo, t) — > are used. Notice that the normalization of 
p is time independent. Indeed, if we write the equation 
in the form dtp = d x J and assume the boundary condi- 
tions ^(±00, t) — ► 0, it can be shown that J_ dx p(x, t) 
is a constant of motion. 

As indicated above, Eq.(Q) can be used to describe a 
large class of anomalous diffusion processes since it con- 
tains, as a particular case, the porous medium equation, 
Levy superdiffusion, as well as a mix of them. For fi = 2 
and v = 1, Eq.(^) recovers the standard Fokker-Planck 
equation with a drift term. The particular case F{x) = 
(no drift), D(x) = constant, p — 2, and arbitrary v has 
been considered by Spohn ]19[ . The case p = 2 has also 
been considered in p^, pl|, |22| with simple drifts. The 
case p < 2 with F(x) = and D — constant has been 
addressed in p3| . Our present analysis involves several 
extensions of these cases, in particular, by employing the 



external force F 



-k\X 



7 . Indeed, a variety 



of systems, which presents anomalous diffusion, can be 
described by short or long-range forces when the previ- 
ous one is employed. This analysis is presented in Sec. 
II for the case p = 2; the solution is obtained in a closed 
form. In Sec. Ill, we also obtain exact results for the 
case p 7^ 2, i.e., for the external drift F — —k\X we ob- 



2 



tain solutions for p G 1Z, and for F = —k\x + /c 7 x|x| 7_1 
we find solutions for p = and /i = 1. Finally, in Sec. 
IV we present our conclusions. 



II. NONLINEAR FOKKER-PLANCK 
EQUATION 

Before starting our discussion, it is interesting to note 
that the results presented in ^FJ, |2l], |2^, |2^, [m| may 
be essentially obtained by using normalized scaled solu- 
tions of the type 



p(x,t) 



1 



*(*)' 



*(*) 



(2) 



in Eq. ([!]). Similar approach has been employed, for in- 
stance, in jH Q. To illustrate this procedure let us 
insert Eq.(@J into Eq.(Q) with p = 2, F{x) — and 
D(x) — T> — constant. Thus, 



d . d 

72— = — 



$ (i) 2 dz 



dz 1 $ (i) 2+,y dz 



(3) 



In the following, we eliminate the explicitly time depen- 
dence on Eq. (ra) choosing 



(4) 



where k is a constant. Thus, Eq. (Q) replaced into Eq. 
(||) gives us 



- k-^-zp(z) = 

ctz az 



^[pwr} • (5) 

From Eq. (^) we obtain 

$(t) = [{l + V )kt] 1/(1+v] , (6) 

where we have adopted the solution which satisfies 
$(0) = for v > — 1. Finally, the solution of Eq. (||) 
leads to 



p(x,t) 



*(*) 



exp„ 



(7) 



with q = 2 — and exp g (x) = [1 + (1 — g)x]— being 
the ^-exponential function. This generalized exponential 
arises within the nonextensive thermostatistical formal- 
ism by optimizing, under appropriate constraints, the en- 
tropic form S q = (l-J dx[p(x,t)]i)/(q - 1) §, |§. The 
constant k can be fixed from the normalization condition 

I-oo dx P( x ' f ) = L 

Let us now go back to Eq. (|]) with F = — k\X + 



A; 7 o;|a;|' 7 1 (fci ^ fc 7 ) and D(x) = T>\x\ 



Firstly, we address our discussion to case p = 2. In this 
case, Eq.(Q) reduces to 



$ (i) d d 



z" & P 



*(*) 



dz 



dz 1 $ (£) 2+1/+s dz 



[p(z)} 1 



z^k^ 



zk\ 



(8) 



when p(x,t), given by Eq. (0), is employed. We are interested in situations for which the scaled solution of the type 
indicated in Eq. (0) is still valid, i.e., we would like to uncouple the t and z dependence in Eq. (|§|) as performed in 
Eq. (H) via Eq. (Ef). This can be accomplished when 7 + 8 + v = 0. If this condition is satisfied, we obtain 



*(t) 



fci 



fc' 



and 



k'^-zp(z) = ^~ 
az dz 



*(t) 2 $(i) 2+l/+S 
d 



vr a -[p(z)r -k^p(z) 



(9) 



(10) 



where A;' is a constant that plays a role analogous to k in Eq. (4), and it is to be determined through the normalization. 
From the solutions of Eq. (|9|) and Eq. (|l0|), which generalize the results of the previous example, we have 



*(*) = 



k' 



:0 



-{x+u+e)kit 



and 



p(x,t) 



exp 



1 



2 + 9 \®(t) 



-1 l/(l+v+9) 



— kry In 



7 1112 -q 



$(t) 



(11) 



(12) 



where ln g a; = (x 1 9 — 1)/(1 — 9) is the q-logarithm func- tion (that is the inverse function of the q-exponential) . 
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In the previous study for the case p = 2, we notice that 
the analytical solution takes fractal and nonlinear 
aspects into account, since there is a spatial dependence 
on the diffusion coefficient {9 ^ 0) and a nonlinear term 
{v 7^ 1). In fact, it contains results obtained in Jl4|, |f^, |2(], 
pij H|, ^5, 27 as particular case. Moreover, it reduces 
to the Rayleigh process j2g] case in the limit q — ► 1 and 
6> = 0. 



III. NONLINEAR FRACTIONAL 
FOKKER-PLANCK EQUATION 

In order to give an extension for the previous discussion 
(fi = 2), we address our investigation concerning Eq. (|l|) 
by considering p 2. Without loss of generality we are 
considering D = 1. We follow the procedure employed in 
fji] , taking the generic property 

^_Q {ax )= a ^G{z) (5 eTZ) (13) 

with z = ax into account. This basic property holds not 
only for the ordinary derivative but also for all fractional 
operators that we are considering here. Thus, substitut- 
ing Eq. (0) into Eq. (Ty) and imposing 



*(*) , fc i 



k' 



$ (i) 2 *(*) $ (t) 



(14) 



in order to uncouple i and z dependence, where k is an 
arbitrary constant. Therefore, we obtain 



*(*) = 



k 2 - — (\ ( " 
k\ 



(15) 



and 

Te — [zp(z) 
az 



d_ 

dz 



d^ 1 

-[p(z)r-k,z- 9 -^+ 2 p(z) 

(16) 



dz^ 



with 7 = —0 — v — p + 2 and k 2 being an arbitrary 
constant. If k 2 = we have $(0) = 0for^ + /i + 6'> 
1, recovering Eq. ([fl]) when /x — > 2. By considering 
&2 = corresponds p(x, 0) oc <5(x). On the other hand, 
employing k 2 ^ leads to p(a;,0) less concentrated than 
5 (x). A more tractable equation to p(z) is obtained after 
an integration of Eq. (pit), i. e., 



k zp(z) = z — — — -j- - A: 7 z 



p(z)+C 



(17) 



where C is another arbitrary constant. In the follow- 
ing, we address our discussion focusing two situations of 
Eq . (|l7|) . The first one is characterized by a linear drift, 
F = —k\x. In the second one, we consider the drift 
F = —k\x + fc 7 a;|a;|' ) ' _1 for particular values of p. 
A. The drift F = -kxx 

Let us start our discussion by considering Eq. (|l7]) with 
k\ ^ and fc 7 = 0. To solve it, we can use the procedure 
employed in |25|| with C = and 



1 +/! 



(18) 



In this case, we have 
.4 



p(x,t) 



$(t) 



(1 + & z )( 1 -^)( 1 +^+ £ ') 



(19) 



with $(£) given by Eq.((T|), z = x/$(t), and 



r(a + i) 



(20) 



Note that Eq. (|19|) reduces that one obtained in |27|] in the 
absence of drift. We notice that the presence of source 
(absorvent) term, apia = const) , in the right hand of 
Eq.(pl) modifies Eq.(p^) as follows: 



$(*) 



(p + 9 



- l)k 



(v-l)a+(p + 6 + v- l)k 



l)ki V 



(i/— l)at 



(21) 



B. The drift F = -kix + k^x^' 1 

Now, we consider other particular cases, namely p = and p = I. Let us start with p = and arbitrary v. The 
corresponding equation is 

^p (x, t) = j*" 9 f [p (y, t)] v dy + [k lX - k 7 x-^ e + 2 ] p(x, t)\ . (22) 
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To solve it let us go back to Eq. ( |T7j ) with p = and C = 0, i.e., 



kzp(z) = z'" / [p (z)] v - k lZ -»- e+2 p(z) 



(23) 



whose solution is given by 



p(z) oc —, = 

k Z l + 6 + fc 7 Z-^+ 2 



I - - 1 7 / dz {k'z 1+e + k^z-^ 2 



1/(1-") 



(24) 



where C is a constant. 

Let us now address our analysis to the case p = 1. It corresponds to investigate the equation 



d_ 

at' 



p(x,t)= x 6 \p(x,t)] y + [kix - k n 



p{x,t) 



(25) 



To obtain the solution of this equation employing the 
ansatz (|2|) we use Eq. (|j~7|) . It follows that 



k zp (z) — z 9 [p(z)Y — k^z 



-v-8+l ~ 



p(z)+C. (26) 



The solution corresponding to C = is 



p (z) oc z 1 " 



(27) 



Let us finally indicate a connection between the results 
obtained here and the solutions that arise from the opti- 
mization of the nonextensive entropy These distribu- 
tions do not coincide for arbitrary value of x. However, 
the comparison between the asymptotic behaviors (large 
\x\) enables us to identify the type of tails. For instance, 
by taking the behavior exhibited in Eq.(n9[) for large \x\ 
with the asymptotic behavior l/|x| 2 /' 9 ~~into account, 
which appears in [|J for the entropic problem, we obtain 



3 + p + e 



1 



(28) 



This relation is the same obtained in |27j in the absence 
of drift. It also recovers that one already established in 
§3| for 9 = 0. 



some space-dependent (power-law) classes of external 
drift and diffusion coefficient. We have shown that it 
admits exact solutions where space scales with a (usually 
simple) function of time (as indicated in Eq.(H)), and 
have fully and explicitly worked out some of those. In 
particular, we first consider the case p = 2 in the pres- 



ence of the drift, F = — kix + k^xlxp^ , with k\, k 1 ^ 0. 
After, we analyse the case p ^ 2 in the presence of the 
same drift for some situations. We have also discussed 
the connection with nonextensive statistics, when appro- 
priated, providing the relation between the entropic in- 
dex q and the exponents appearing in the Fokkcr-Planck 
equation. Thus, by extending here results such as those 



discussed in 23 24 



25J, |27|, we hope to make possible 
further applications to physical systems exhibiting suit- 
able mix of different forms (nonlinear v ^ 1, fractional 
derivative p ^ 2 and fractal 6 ^ 0) of anomalous diffu- 
sion. 
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